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Abstract 

We present an updated analysis of the current solar neutrino data in terms of three-flavor os- 
cillations, including the additional constraints coming from terrestrial neutrino oscillation searches 
at the CHOOZ (reactor), Super-Kamiokande (atmospheric), and KEK-to-Kamioka (accelerator) 
experiments. The best fit is reached for the subcase of two-family mixing, and the additional ad- 
mixture with the third neutrino is severely limited. We discuss the relevant features of the globally 
allowed regions in the oscillation parameter space, as well as their impact on the amplitude of 
possible CP-violation effects at future accelerator experiments and on the reconstruction accuracy 
of the mass-mixing oscillation parameters at the KamLAND reactor experiment. 

PACS numbers: 26.65.+t, 13.15.+g, 14.60.Pq, 91.35.-x 
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I. INTRODUCTION 



In this paper we present an updated three-flavor oscillation analysis of the solar neutrino 
data coming from the Homestake 0, SAGE @, GALLEX/GNO || §, Super-Kamiokande 
(SK) ||, and Sudbury Neutrino Observatory (SNO) || 0] experiments. We include ad- 
ditional constraints coming from terrestrial neutrino oscillation searches performed at the 
CHOOZ reactors M, at the SK atmospheric v experiment [fj, and at the KEK-to-Kamioka 



(K2K) long baseline accelerator experiment [10|. Implications for upcoming or future ex- 



periments are also discussed. This work extends a previous solar 2v analysis JTT], to which 
we refer the reader for technical details. 

The plan of our paper is as follows. In Sec. II and III we briefly review the theoretical 
and experimental input, respectively. In Sec. IV we present the results of our 3z/ oscillation 
analysis. In Sec. V we discuss some implications for CP violation searches at future accel- 
erator experiments. In Sec. VI we analyze the accuracy of parameter reconstruction in the 
Kamioka Liquid scintillator Anti Neutrino Detector (KamLAND) experiment at reactors 



12 . We draw our conclusions in Sec. VII. 



II. THEORETICAL 3i/ FRAMEWORK AND APPROXIMATIONS 

We consider standard 3u oscillations among flavor eigenstates v a = (u e , u T ) and mass 
eigenstates z/j = u 2 , v 3 ) with a squared v mass spectrum defined as 

(Sin? Std? \ 
2~' + — ±Am2 J ' W 

up to an irrelevant overall constant. The two squared mass gaps 5m 2 and Am 2 (both > 0) 
are usually referred to as "solar" and "atmospheric" squared mass gaps, respectively. The 
cases +Am 2 and — Am 2 identify the so-called direct and inverted spectrum hierarchies ||13|| . 1 
The mixing matrix XJ ai is parametrized according to the usual convention 



as 



u = u(e 12 ,e 13 ,e 23 ) , (2) 

where we have dropped the (currently unobservable) CP violation phase. The mixing angles 
span two octants, dij G [0, 7r/2], and are often parametrized in terms of either tan 2 (in 
logarithmic scale) or sin 2 (in linear scale) (15], [16| . We use the latter representation in 
this work. 2 

As it is well known, solar and terrestrial neutrino data consistently favor the so-called 
hierarchical hypothesis 

5m 2 < Am 2 , (3) 

which can be used to simplify the calculations JTJ]]. Even if some regions of the explored 
parameter space do not fulfill the above hypothesis, the a posteriori likelihood of combined 



1 We assume direct hierarchy, unless otherwise stated. As discussed later, the case of inverse hierarchy does 
not appreciably change the results of our analysis. 

2 The alternative representation in terms of log tan 2 9ij is more useful when the allowed regions span several 
decades in mixing angles — a situation no longer realized in the current neutrino phenomenology 
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solutions in such regions appears to be so low that, at present, the zeroth-order approxima- 
tion (one- mass-scale dominance [15, 16, [T^] ) is practically justified in the current analysis of 



both solar and atmospheric neutrino data. 

In particular, solar neutrino oscillations are described with very good accuracy in terms of 
the 3u parameter subset (5m 2 , # 12 , #13) ]TE|, [IB], [HJ. Corrections due to violations of Eq. ([|) 
have been shown to be typically small (see [T3|] and references therein), and can be safely 



neglected at present. Therefore, in the statistical analysis of the solar v data, the x 2 function 
takes the form: 

xLar = xWA2,M , (4) 

with no dependence on Am 2 and thus on the mass spectrum hierarchy (direct or inverse). 

Concerning atmospheric neutrinos, it has been shown in several analyses that the usual 
bounds on the dominant — > v T oscillation parameters (Am 2 , #23) axe rather robust under 



perturbations beyond the one-mass-scale dominance, induced by either 13 > [g, |20], |21], |22 
or by 5m 2 > or both [^9, |4], |25|, |26| (within reactor bounds). 3 



Moreover, since the 623 value is irrelevant for the analysis of both solar and reactor 
neutrino data, we only need to know the atmospheric Am 2 scale in our analysis, as derived 
by a combination of the (preliminary) SK and K2K bounds || |llj for unconstrained 623 (see 
the next section). Summarizing, we use the SK (atmospheric) and K2K bounds to derive 
the marginal likelihood for Am 2 , namely, 

XSK+K2K = XsK+K2K(^ m2 ) ■ (5) 

Such likelihood is practically independent of all the other 3u parameters (and on the ±Am 2 
cases), as far as the hierarchical hypothesis holds phenomenologically. 

In the analysis of the CHOOZ reactor constraints, however, the hierarchical approxima- 
tion may not be accurate enough, and the full 3z/ oscillation probability (as reported, e.g., 
in [[n|) must be used in global analyses |23|, ^Sj, |2D|, |3D], providing a x 2 functional 



dependence of the kind 

Xchooz = XchoozO™ 2 > ±Am 2 , 12 , 9 n ) . (6) 

The reason is that the CHOOZ bounds on V e disappearance can be saturated by taking 
either #13 > or 5m 2 > or both [ER 1261 EM |32| , thus providing an anticorrelation between 



the CHOOZ upper limits on 5m and #13 that must be properly taken into account. In other 
words, the higher 5m 2 , the lower 813, in both cases of normal and inverted hierarchy [^, |29 . 
We will discuss some consequences of such anticorrelation in Sec. V. 

Before performing the final combination of SK+K2K+CHOOZ with solar neutrinos, it 
is useful to project away the functional dependence of the likelihood on Am 2 by defining a 
"terrestrial" x 2 as 

X t 2 err (W, 9 12 , Ota) = mm f Xsk+K2K + Xchooz) , (?) 

which contains only a residual dependence on the hierarchy (direct or inverse) through 
Xchooz- Therefore, in principle, the global combination 

Xglobal^ 7 ™ 2 ' 012, 013) = Xsolar + Xterr ( 8 ) 

also depends on the hierarchy, although, unfortunately, in a very weak way at present. 



3 Moreover, the estimates of the leading parameters (Am 2 , 623) are robust under perturbations induced by 
new neutrino states or interactions, see p7| and references therein. 
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III. EXPERIMENTAL INPUT 



In this section we briefly review the experimental input coming from solar and terrestrial 
(atmospheric, reactor, and accelerator) experiments. 

For the solar neutrino analysis we use the same standard solar model input and neutrino 
data (81 observables) as in [jO}] , but with an updated winter-summer rate difference from 
the combination of GALLEX/GNO § 4 and SAGE § measurements, 

R w - R s ~ -5 ± 7 SNU (GALLEX/GNO + SAGE) . (9) 

Concerning reactor data, we take the CHOOZ absolute spectrum from ||. The 3v analysis 
includes the 7+7 bin spectra as in fTB |. 



Concerning atmospheric neutrinos, the SK collaboration has recently presented updated 
(92 kton yr) bounds on the leading — > v T atmospheric neutrino oscillation parameters 
(Am 2 , sin 2 26* 23 ) ||. Since the correlation between the two parameters appears to be negli- 
gible, the bounds on Am 2 for unconstrained (i.e., projected) values of sin 2 2#23 are basically 
obtained by fixing sin 2 26*23 at its the best fit value (unity). From a graphical reduction of 
the bounds shown in 0, we derive then an approximate Ax| K function in terms of Am 2 . 

The K2K Collaboration has also presented preliminary bounds on the same mass-mixing 
parameters (Am 2 , sin 2 2# 23 ) UTO]] , again with apparently negligible correlation around the 
best fit point (as far as sin 2 2# 23 is nearly maximal). The bounds on sin 2 2# 23 are relatively 
weak, but those on Am 2 are already competitive with those placed by SK and must be taken 
into account. We simply do so by a graphical reduction of the Ax|- 2K function presented in 
|T0| for maximal mixing, to be combined (summed) with the previous Axg K . 

It turns out that the function Axg K+K2K obtained in this way is nearly parabolic in the 
linear variable Am 2 , although the separate Ax| K and components are not parabolic. 

The SK+K2K x 2 minimum is at 2.7 x 1CT 3 eV 2 , with a ±la error (Ax1k+K2K = 1 f° r 
Ndf = 1) equal to ~ 0.4 x 1CT 3 eV 2 . The estimated range for Am 2 , as discussed in the 
previous section, is known to be very stable under small perturbations induced by nonzero 
values of #13 or 5m 2 . Therefore, we use the following SK+K2K combination, 

Am 2 ~ (2.7 ± 0.4) x 10~ 3 eV 2 (10) 

with approximately linear, symmetrical and gaussian errors, for unconstrained values of the 
other 3v parameters. Notice that values of Am 2 below 1.5 x 10~ 3 eV 2 are excluded at 



3a, corroborating the trend of a previous SK+K2K combination |33| , and reinforcing the 
validity of Eq. ([3|). Of course, it will be desirable to confirm this nice result by undertaking 
a thorough and ab initio analysis and combination of the K2K and SK data. We plan to do 
so when such preliminary data will be published and described in more detail. 5 



IV. RESULTS 

In this section we describe the main results of our 3v solar neutrino oscillation analysis 
with additional terrestrial constraints. We start from the subcase of 1v oscillations. 



4 Corrected to remove eccentricity effects. 

5 In particular, the K2K spectral shape analysis leading to the results in fiofl appears to involve much more 
(and currently unpublished) information than it was required in the simpler analysis of the total rate only 
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TABLE I: Absolute and local best fits for the 2v analysis reported in Fig. 1 (without CHOOZ). 
The value Xmin = 72.3, reached in the LMA solution, corresponds to a good overall fit to 81 solar 
neutrino data (minus 2 free parameters). 



Solution 


5m 2 (eV 2 ) 


sin 2 6 12 


Ax 2 


LMA 


5.5 x 1(T 5 


0.30 




QVO 


6.5 x 1(T 10 


0.55 


7.5 


LOW 


7.3 x 1(T 8 


0.41 


10.5 



A. 2v oscillations 

In the subcase of pure 2v oscillations (#13 = 0), the information on Am 2 coming from 
SK+K2K is completely decoupled within our approximations, and the solar+CHOOZ v 
parameter space reduces to (5m 2 , #12). 

Figure 1 shows the constraints coming from the global solar neutrino analysis in the 
(5m 2 , sin 2 #12) plane, with and without CHOOZ. In the case of solar v data only, the solutions 



are almost coincident with those presented in Fig. 1 of |TT[ , modulo the small changes induced 
by the updated input from Eq. @ and by the different abscissa. For completeness, Table I 
gives the local \ 2 minima for the so-called large mixing angle (LMA), quasivacuum oscillation 
(QVO) and low-£m 2 (LOW) solutions. These minima are not altered by the inclusion of 
the CHOOZ data, whose only effect is to slightly strengthen the upper bounds on 5m 2 , as 
shown in Figure 1. The features of such 2v allowed regions have been discussed in a number 
of papers B U\, B3, K5l EB, [37], BR |3!| and are not repeated here. 



B. ?>v oscillations 



Figure 2 shows sections of the volume allowed by solar neutrino data only in the 3z/ 
parameter space (5m 2 , sin 2 # 12 , sin 2 13 ), for four representative values of sin 2 9% 3 (equal to 
0, 0.02, 0.04, and 0.06). The best fit is reached in the same LMA point as for the 2v case 
(sin 2 #13 = 0). The constraints are derived through Eq. (f|) and A^df = 3. It can be seen that, 
for increasing sin 2 #13, the upper bounds on 5m 2 in the LMA become slightly weaker, and 
the LOW and QVO solutions become less unlikely. Maximal (v\, 2/2) mixing (sin 2 #12 = 1/2) 
is also less disfavored for increasing sin 2 6*13. This trend is qualitatively consistent with the 
recent results in |35| (obtained for an almost equivalent solar v data set), as well as with 
previous "iv analyses of solar neutrino data |L8], ^T], BJ], |IJ. 

Looser constraints on either 5m 2 or maximal mixing can be a desirable property from a 
theoretical or experimental point of view. Relatively large values of 5m 2 are, e.g., required 
for the possible detection of leptonic CP violation in future (very) long baseline experiments 
(see, e.g., [P]]). Nearly maximal mixing can allow certain cancellations in 0z/2/3 decay or in v 
model building (see, e.g., [22]]). Unfortunately, such "desirable" features in Fig. 2 are severely 
limited by the inclusion of current terrestrial data through Eqs. (|7|) and ([5]). This is shown 
in Fig. 3, where the addition of the terrestrial v constraints actually strengthens the upper 
bounds on 5m 2 for increasing sin 2 # 13 , since the CHOOZ bounds are more quickly saturated. 
Notice also that, in comparison with earlier analyses, the current K2K+SK bounds on Am 2 
[Eq. (|10D] are now sufficiently strong to prevent the occurrence of relatively low values of 
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Am 2 (< 1.5 x 10~ 3 eV 2 ), where the CHOOZ bounds would become significantly weaker. 
Therefore, the current SK+K2K data enhance the impact of terrestrial constraints on the 
solar 3u parameters (5m 2 , sin 2 #i 2 , sin 2 6*13). 

Figure 4 provides another way of looking at the results of the solar+terrestrial analysis. 
In this figure, the Ax 2 function is projected onto one parameter at a time, so that the 
n-sigma bounds on each of them (the others being unconstrained) are simply given by a 
Ax 2 = n 2 cut. In the left panel it can be seen that the x 2 difference between the LMA 
and the LOW or QVO solutions is slightly lowered in the Zv fit, as compared with the 2v 
case reported in Table I. However, such improvement for the LOW and QVO solutions for 
nonzero #13 is currently marginal (less than one unit in Ax 2 ). Basically all data converge 
towards zero or low values of sin 2 #13, whose global upper bounds are reported in the right 
panel of Fig. 4. In particular, one gets 

sin 2 #i3 < 0.05 (3cr, N DF = 1) , (11) 

for unconstrained values of the other 3z/ parameters. 

In conclusion, genuine 3u mixing in solar neutrinos is severely constrained by terrestrial 
data. For the most favored solution, increasing #13 corresponds to lowering the upper bound 
on 8m 2 . For the less favored LOW and QVO solutions, the improvement of the fit for slightly 
nonzero values of 9 13 (see also [p8f] ) is not statistically significant. Figures 3 and 4 quantify 
such effects. 

A final remark is in order. Figures 3 and 4 refer to the case of direct neutrino spectrum 
hierarchy. We have verified that the changes for inverse hierarchy (not shown) are negligible. 
Indeed, the phenomenological difference between the two hierarchies, which vanishes for 
#13 —> 0, is still too small to emerge for the nonzero values of #13 allowed by Eq. (|Tl|). 



V. IMPLICATIONS FOR FUTURE CP VIOLATION SEARCHES 



The anticorrelation between the upper bounds on 5m 2 and sin 2 # 13 shown in Fig. 3 is not 
a particularly desirable feature for future CP violation searches at accelerators, where the 
difference between the v and V oscillation probabilities would be enhanced for both 5m 2 and 
#13 close to their upper limits (see, e.g., the reviews in [fTI|j). 

In order to study in more detail the maximum CP violation effects compatible with 
current bounds, we think it useful to focus on the following benchmark quantity, which 
typically appears as an overall prefactor modulating CP-violating amplitudes, and which 
has the useful property to be a vacuum-matter ( VM) invariant [^3| : 



'VM 



(ml - m 2 2 ) (m 2 2 - m 2 )(m 2 - m\) sin2#i 2 sin2#i 3 cos# 



13 



5m 



4(Am 2 ) 2 - (dm 2 ) 2 



sin #12 cos #1 2 sin 3 cos 9 



'12 



13 1 



'13 



(12) 
(13) 



More precisely, in order to eliminate the dependence on Am 2 and to simplify the discussion, 
we integrate /vm over the gaussian distribution in Am 2 defined in Eq. (|II]), obtaining the 
expectation value 



■ VM 



5m 



A(Am 2 Q ) 2 + 4(<7 A ) - (5m ) sin 6 12 cos 9 12 sin #13 cos 2 #13 



(14) 



where Am 2 = 2.7 x 10" 3 eV 2 and a A = 0.4 x 10" 3 eV 2 [Eq. (0)]. The quantity I yu is then 
a function of the same 3u parameters (5m 2 , #12, #13) used for the previous solar+terrestrial 
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analysis. Notice that 7vm — > when any of these 3z/ parameters tends to zero, as intuitively 
expected for a quantity associated to CP violation flS ]. 

Figure 5 shows isolines of Jvm in the plane (5m 2 , sin 2 # 12 ), for the same representative 
values of sin 2 # 13 as in Figs. 2 and 3. The 5m 2 scale is restricted to the range relevant for 
the LMA solution, whose C.L. contours are superposed in each panel (as taken from Fig. 3). 
In the upper left panel (sin 2 #13 = 0), the vacuum- matter invariant Jvm is identically zero. 
In other panels, it can be seen that Jvm achieves its maximum values, within the LMA 
solution, for "intermediate" values of sin 2 #13 (~ 0.02-0.04). Indeed, for increasing values of 
sin 2 the LMA upper bound on Sm 2 decreases, and Jvm is maximized when a compromise 
is reached. We conclude that the maximum allowed amplitude of possible CP violation 
effects (in future accelerator experiments) occurs somewhat below the upper bound placed 
by terrestrial experiments on sin 2 6^3, due to the phenomenological anticorrelation with the 
5m 2 upper bound. 



VI. IMPLICATIONS FOR KAMLAND 



The KamLAND experiment is currently detecting reactor V e events, whose energy spec- 
trum will provide soon a clear terrestrial test of the solar v LMA solution [12]. If the LMA 
solution is disconfirmed, the BOREXINO experiment will be able to test the remaining 
LOW and QVO solutions by detecting time variations of the event rate [53]. Assuming that 
the LMA solution is confirmed, various studies have shown how well the KamLAND exper- 
iment can reconstruct the oscillation parameters for a few selected points in the (Sm 2 , 612) 
parameter space |T2L H3, [19, E7), |48], [41], ^0], |5T|. In this section we extend these studies by 
performing a continuous scan in the 5m 2 parameter, in the region of interest for the LMA 
solution. 



In the analysis, we use standard inputs for the reactor average fuel composition |p2 
energy spectra 0] and cross section ||54|| . We take from |55| the reactor power and locations. 
The visible energy window is chosen to be E e [1.22, 7.22] MeV (divided in 12 bins) as in 
fH, H, with energy resolution g(E)/VE = 5% (E in MeV) @. We assume a "KamLAND 
year" corresponding to 550 events/yr for no oscillation |12|, and consider two representative 



exposure periods of 0.5 and 3 years, the shortest period being a tentative estimate of the 
detector live-time that might be used for the first official release of the KamLAND data. 

Concerning the uncertainties, in addition to the statistical fluctuations in each bin, we 
include a ~ 4% overall normalization error |I2| (due to uncertainties in fiducial volume, v 
flux, and cross section), and a 2% energy scale uncertainty fL2| . The normalization and 
energy scale systematics (fully correlated in each bin) mainly worsen the reconstruction 
accuracy of sin 2 #12 and 5m 2 , respectively. Background reduction and subtraction are still in 
progress at KamLAND ||12|| , and we do not include the (unknown) associated errors in the 
analysis. More realistic analyses, including backgrounds and other systematics (e.g., fuel 
composition uncertainties will be possible after release of the real KamLAND data 

and of related information. 

Given the above input, for fixed true values of the parameters (5m 2 , sin 2 #12) we first 
determine the Ax 2 < 4 region including the reconstructed oscillation parameters, and then 
project this region either on 5m 2 or on sin 2 ^, in order to obtain the ±2<r error bands 
on each of the two oscillation parameters separately. For simplicity, this procedure is done 
for fixed #13 = 0. When real KamLAND data will be available, a more complete analysis 
procedure (not simulated here) should include a variable #13 parameter, to be constrained 



7 



by world neutrino data (as previously done for the currently available data in Figs 3 and 4). 
Small variations of #13 would then basically act as an additional KamLAND normalization 
uncertainty |49| . 

Figure 6 shows the ±2<r range of the reconstructed 5m 2 , as a function of the true value of 
5m 2 . The true value of sin 2 # 12 is fixed at 0.3 (LMA best-fit), while the reconstructed sin 2 # 12 
is projected away. It can be seen that, in the range most relevant for the LMA solution 
(5m 2 ~ few x 10~ 5 eV 2 ) the reconstruction accuracy depends only weakly on 5m 2 . It become 
rapidly worse, however, outside this range. E.g., for low values of 5m 2 (approaching the no- 
oscillation limit in KamLAND), the reconstructed value of 5m 2 can be considerably higher 
than the true one, the variation being compensated by a lower value of the reconstructed 
sin 2 #12. For high values of 5m 2 (say, ^ 1CT 4 eV 2 ), degenerate solutions begin to appear, 
eventually leading to a tower of (merging) multiple ranges for the reconstructed 5m 2 . For a 
half-year exposure, an isolated, low-<5m 2 reconstructed range also appears, since the spectrum 
suppression due to true values 5m 2 > 10~ 4 eV 2 can be confused with a similar suppression 
induced by lower values of both 5m 2 and sin 2 # 12 . Such ambiguities can be cleared up, 
e.g., after 3 years of data taking, provided that 5m 2 £ 2 x 10~ 4 eV 2 . Above this range, 
KamLAND cannot resolve the pattern of fast oscillations in the I7 e energy spectrum, which 



could be disentangled only by reducing the typical baseline [30, 50, 56 



Figure 7 zooms in on true 5m 2 values below 2 x 10 -4 eV 2 , for three representative (true) 
values of sin 2 #12. For oscillation parameters close to the LMA ones, the accuracy in the 
5m 2 reconstruction can already be as good as ~ ±10% at 2a after half-year exposure, and 
can improve by a factor of ~ 2 after three years. The 5m 2 reconstruction accuracy improves 
(worsens) slightly for increasing (decreasing) values of sin 2 6*i 2 , since this parameter governs 
the amplitude of the oscillation pattern, from which 5m 2 is inferred. 

Finally, Fig. 8 shows the reconstructed range of sin 2 #12 (first octant only), for the same 
choice of true oscillation parameters as in Fig. 7. From this figure it appears that the most 
accurate reconstruction of sin 2 #12 occurs for 5m 2 values slightly below the current LMA 
best-fit point. The reason is that, for such relatively low values, the overall rate suppression 
in KamLAND is maximized, making it easier to put bounds on the oscillation amplitude 
and thus on sin 2 # 12 . The sin 2 # 12 reconstruction accuracy tends to become constant for high 
values of 5m 2 (regime of fast oscillations). In the lower panel, notice the appearance of a 
degenerate reconstructed range for the case of 0.5 year exposure, which has the same origin 
as the corresponding one discussed in Fig. 6. 

In conclusion, KamLAND is in an extremely favorable situation to fix the most likely 
values of 5m 2 in the LMA region within a few % error. The achievable accuracy in the 
sin 2 612 reconstruction is instead both lower and more uncertain, and will depend very much 
on how the real KamLAND data will look like. 



VII. CONCLUSIONS 

We have studied the 3z/ perturbations to the usual 2v solutions to the solar neutrino prob- 
lem, including — besides the full set of current solar data — the terrestrial neutrino constraints 
coming from the CHOOZ (reactor), SK (atmospheric) and K2K (accelerator) experiments. 
The global best fit is reached in the so-called LMA solution, for the subcase of 2u mixing 
(#13 = 0). The likelihood of the (less favored) LOW and QVO solutions does not improve 
significantly for nonzero #13, within the severe constraints placed by terrestrial neutrino ex- 
periments on such mixing angle. The anticorrelation between the upper bounds on 5m 2 
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and 6 is has been discussed, together with its implications for the maximum allowed ampli- 
tude of possible CP violation effects (expressed in terms of a representative vacuum-matter 
invariant). We have also discussed the expected accuracy of the mass-mixing parameter 
reconstruction in the (currently running) KamLAND experiments, described through a con- 
tinuous scan of the Sm 2 parameter within the LMA solution. 
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2v oscillations! all solar data w/o CHOOZ 
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FIG. 1: Two-flavor global analysis of solar neutrino oscillations in the (5m 2 , sin 2 612) parameter 
space, with and without the additional constraints placed by the CHOOZ reactor experiment. The 
inclusion of CHOOZ leads to slightly more restrictive upper bounds on dm 2 . 
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FIG. 2: Three-flavor global analysis of solar neutrino oscillations in the (6m 2 , sin 2 #12, sin 2 #13) 
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parameter space, shown through four sections at fixed values of s 2 3 = sin 2 0\ 
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3v oscillations 1 Solar + Terrestrial data 
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FIG. 3: As in Fig. but including the constraints from terrestrial neutrino oscillation searches 
at CHOOZ, SK, and K2K. 
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FIG. 4: Projections of the global (solar +terrestrial) Ax 2 function onto each of the 
(5m 2 , sin 2 #12, sin 2 #13) parameters. The n-sigma bounds on each parameter (the others being 
unconstrained) correspond to Ax 2 = n 2 . 
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FIG. 5: Isolines of the vacuum-matter invariant Ivm, m units of 10 eV . Superposed are the 
contours of the large mixing angle solution. The invariant Ivm appears as a typical prefactor in 
CP violating amplitudes expected at future accelerator searches; see the text for details. 
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KamLAND reconstruction of <5m 2 at ±2a 
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FIG. 6: Simulation of KamLAND results. Reconstructed range of 5m 2 at ±2cr, as a function of 
the true value of 5m 2 in abscissa. The true value of sin 2 9\2 is fixed at 0.3. The curves refer to 
detector exposures of 0.5 and 3 years. 
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KamLAND reconstruction of <5m 2 at ±2a 
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FIG. 7: Simulation of KamLAND results. Reconstructed range of 5m? at ±2<r (normalized to the 
true value of 5m? in abscissa), for three representative values of sin 2 #12, and for exposures of 0.5 
and 3 years. 
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KomLAND reconstruction of sin 2/ # 12 at ±2o~ 
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FIG. 8: Simulation of KamLAND results. Reconstructed range of sin 2 #12 at ±2cr as a function 
of 5m 2 , for three representative (true) values of sin 2 612, and for exposures of 0.5 and 3 years. 
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